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Abstract 

We analyze a one-dimensional quantum model with off-diagonal disorder, consisting of a sequence of 
potential energy barriers whose width is a random variable either uniformly or "half-normally" distributed, 
subjected to an external electric field. We shed light on how the microscopic disorder affects the value of 
the transmission coefficient, and on the structure of the fluctuations around the solutions corresponding to 
the regular lattice configuration. We also characterize the asymptotic limit obtained by letting the number 
of barriers diverge, while their total width is kept constant. Thus, we explain the novelty of our method 
with respect to the standard thermodynamic limit discussed in the literature, and also evidence the onset of 
a large deviations principle for the transmission coefficient. 
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I. INTRODUCTION 



Nonequilibrium thermodynamics is based on the notion of space and time scales separation and 



on the assumption of local equilibrium The theory of large deviations, in particular, helped 



to understand and interpret the role of the fluctuations in nonequilibrium systems HVNlOll. On the 
other hand, recent technological advances on the nanoscale science and technology demand an 
extension of the theoretical apparatus and foster a statistical mechanical approach to systems of 
relatively small numbers of degrees of freedom. In such systems the microscopic, mesoscopic 
and macroscopic scales can not be sharply separated, and the physical properties of microscopic 
devices widely fluctuate with respect to their mean values, violating the standard thermodynamic 
laws which describe macroscopic fields. In this work we face these issues by considering a variant 



of the original Anderson model, which is the prototype of a disordered solid DUD- In particular, we 
investigate the role of the microscopic disorder on the transmission coefficient of one-dimensional 
systems consisting of a sequence of N barriers, with random widths, and — 1 wells, under the 
constraint that the sum of the barrier widths and the total length of the system are fixed and do not 
change with A^. We then introduce a classical thermostat at given temperature T and an external 
electric potential V^ — Vr. Furthermore, we do not introduce simplifying assurnptions such as the 



ization 



"tight-binding" approximation introduced by Anderson in his pioneering paper [|12| 1 on loca 
effects in disordered solids. Therefore, our model enjoys a purely off-diagonal disorder [[IS 
which concems only the tunneling couplings among the wells, leaving unaffected the energies 
of the bound states within the wells. This is not the case of the original tight-binding model, 
whose random fluctuations only concern the energy of a bound state. In turn, while in Anderson's 
model increasing the number of barriers corresponds to taking the large system limit, in our case 
it corresponds to distribute more finely the same amount of insulating material within the fixed 
length of the system. 

The introduction of an external field allows to extend to "nonequilibrium" the results previously 



obtained in the analysis of the model treated in [|l6|] . which are recovered, as shown below, in the 



limit of vanishing external fields. We, thus, investigated the effect induced by this kind of disorder 
at the mesoscopic scale on the transmission coefficient and we shed light on the structure of its 
fluctuations. Our results can be summarized as follows: 

• There are no localization effects for the equilibrium distribution of energies at temperature 
T = 300i^': positive currents persist even in the large A^ limit. 



Furstenberg type theorems flllH do not apply. The reason is that the product of the random 
matrices yielding the transmission coefficient for a given choice of barriers changes, in 
order to preserve the length of the system and the sum of the barrier widths, when the + 1- 
th barrier is introduced. 

The value of the transmission coefficient, averaged over an ensemble of disordered configu- 
rations, is close, for large N, to the value corresponding to the ordered sequence of equally 
spaced barriers and wells, which is bounded away from zero. 

There is a scale for A^, above which the (always positive) transmission coefficient does not 
depend on the specific realization of the disorder, but still depends on A^, and there is another 
scale above which even the dependence on A^ is eliminated. We call "mesoscopic" the first, 
and "macroscopic" the latter scale, since it represents macroscopic nanostructured materials. 
This means that all realizations of the disorder become equivalent in the N oo limit. 

At room temperature, the probability distribution function (PDF) of the time independent 
transmission coefficients of the different realizations of the system satisfies a principle of 
large deviations. Furthermore, the peak of this PDF corresponds to the transmission coeffi- 
cient of the regular realizations. 

Our N oo limit, representing a macroscopic object at given temperature, which is mi- 
croscopically randomly structured, leads to radically different results from the usual macro- 
scopic limits. In particular, it leads, in certain systems, to the experimentally verifiable 
lack of localization. This is relevant in situations complementary to those described by the 
standard theories. 



II. THE MODEL 

Our one-dimensional model of a macroscopic semiconductor device consists of an array of A^ 
potential barriers and A^ — 1 conducting regions (wells), in contact with one electrode which acts 
as an extemal thermostat at temperature T = 300K. The particles leaving this thermostat are 
subjected to an extemal electric field F, cf. Fig. [U The barriers have a constant height V{x) = V 
while their width is either uniformly or "half-normally" randomly distributed. For any A^, the 
widths of the conducting regions take a constant value 5^- We denote by L the fixed total length 
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of the sample, by the fixed sum of the widths of all the barriers (i.e. the total length of the 
insulating region), and by /3 the fixed ratio between insulating and conducting lengths, so that 

L = {l + f5){N -1)5n (1) 
holds. To compute the current, we study the steady state Schrodinger Equation (SE): 

—i)ix) + iV-eFx)i)ix)=Ei)ix), xG[0,L] (2) 

2m dx'^ 

where m is the mass of the particle, e the electronic charge and F the magnitude of the external 
electric field which takes the values F''^ inside the barriers and in the conducting regions. Due 
to the electric field, the potential energy decreases monotonically from Vi, on the left boundary, 
down to Vr, on the right boundary, with a slope given by, respectively, —eFis within the barriers 
and — eFc in the conducting regions. Let us also introduce the parameter r = Fig/F^ < oo, which 
allows to consider the presence of a nonvanishing electric field even within the wells. Therefore, 
the energy of the electric field acting on the system, denoted by Ey, amounts to 

Ey = e{Ve - Vr) = eF^rP + 1) (3) 

The boundary conditions prescribe Aq > for the amplitude of the plane wave entering from 
the left boundary and A^^^i = (no wave enters or is reflected from the right boundary). The 
barriers are delimited by a set of 2N points, denoted by xq = 0, X2n-i = L in Fig. [H hereafter 
called nodes of discontinuity of the potential. The left boundary consists of a classical thermostat 
at temperature T, from which particles emerge at node in xq as plane waves, with energies dis- 
tributed according to the Maxwell-Boltzmann distribution. Differently, no particles come from the 
electrode on the right. 

Thus, denoting by Ue the i-th region, for i E {0, 2, 2A^}, the solutions of eq. ^ take the 
form: 

Mx) = A2eAii{x) + A2i+iBit{x), ifxeUi (4) 

where Aii{x) and Bii(x) denote the Airy functions. In each of the conducting regions, one may 
define the steady state currents as: 



Je{x) 



2mi 



dx J dx 



j^iA^e) - j f {A^e+i) , (5) 




FIG. 1: (Color online) 1-D multiple-barrier system, consisting in a sequence of 2N + 1 regions: barriers, 
whose width is uniformly random distributed, (A^ — 1) conducting regions of constant width 5n, and the 
two boundary regions (shadowed areas), on the left and on the right, representing the electrodes. The latter 
are modelized as a classical thermostat at temperature T = SOOJC and electric potential Vi (electrode on the 
left) and Vr (electrode on the right). The mean energy of the plane waves entering from the left boundary is 
given by ^ksT. The picture illustrates the case characterized by vanishing electric field in the conducting 
regions. 

where the * denotes complex conjugation, f/'{A2e) = hk/m\A2e\'^ denotes the current transmitted 
from the (£— l)-th barrier on the left (or, for i = 0, from the thermostat located at the left boundary) 
and jf" {A2i+i) = /ifc/m|y42£+ip denotes the current reflected from the (i + l)-th barrier, cf. Fig. 
[TJ Then, the application of the BenDaniel-Duke boundary conditions [17], which require the 
continuity of -ipeix) and -^^iix) at the nodes, results in the constancy of the value ji{x) across the 
wells and entails je{x) = j2N{x), for every even Equation (|5]), together with Eq. (HJ, leads to 
the following definition of the transmission coefficient J across the system: 

which depends on the several parameters of the model, such as the number of barriers A^, the 
energy of the applied electric field, the ratio r and the disordered configuration of the sequence 
of barriers. In Sec. |V]we focus, in particular, on the structure of the fluctuations of J as a function 
of the realization of the disorder. We also discuss how the magnitude of these fluctuations depends 
on N, for a given distribution of noise realizations, by exploring a large range of scales: from 
the microscopic one, where = 0(1), up to the macroscopic one, with A^ ^ 1. In order to 
numerically compute the coefficient J as a function of the various parameters of the model, it 
proves convenient to rescale Eq. (|2]l with respect to characteristic quantities, in order to rewrite 
it in a dimensionless form. For this purpose, let us introduce x = Lx, with L given by (O, 
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E = EEt, ^{x) = V = VEt, = F^ET/{eL), with Et = keT (i.e. twice 

the mean kinetic energy of the plane waves entering the bulk from the left side). Moreover, by 
introducing the scalar parameter a = h?/ {2mL'^ET), one obtains the following expression for the 
dimensionless wavevectors: k = ^JoL\fE and z = ^Ja\/'V — E. In the sequel we will refer to the 
dimensionless quantities and, to this aim, we may omit the tilde symbols, for sake of simplicity. 
The dimensionless version of eq. (|2]l, then, attains the form: 

-a—^{x)^iy{x)-eFx)^{x) = E^{x), a; G [0, 1] (7) 

and it represents the SE which will be solved numerically with the aforementioned conditions at 
the nodes. 



ni. THE TRANSFER MATRIX TECHNIQUE 



Let us now describe our method of solution of the SE, eq. dV]), which follows Refs. rtllL llTI] 
and is referred to as the Transfer Matrix (TM) technique. Using eqs. (01), the BenDaniel-Duke 
boundary conditions on the generic £-th node, with i G {0, 1, 2N — 1}, cf. Fig. [H read as: 

where = Yll=i + (^f if ^ is even and xe{x) = Yjf=i^^'^ ^« + if ^ i^ odd and where Aj 
denotes the random width of the i-th barrier. Thus, in matrix form, eqs. dH) takes the form: 





(9) 
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where the 2x2 matrices of coefficients 'WL2e{xe) and M2^+i(x^) read 



M„(.,) = I ^"}''^ ^':<^'> I (10) 

and 

M2e+i{xe) =1 ,1 (11) 

\Ai^^-^{xe) Bi^^-^{xi) 

Assuming that the amplitude Aq of the incoming wave is known and that A^n^i = because 
there is no reflection from the right boundary in the 2A^-th region, Fig. \T\ then the linear system 
Q corresponds to a set of AN equations in the AN variables {Ai, ...,A^n}. Skipping the spatial 
dependence of the matrices to keep our notation simple, eq. dH) implies 

11 1 / A^N \ I ^AN \ 

Mq^ ■ Ml ■ ^ ■ Ms ■ ■■■ ■ M,^_, ■ M4Ar-i^ - I ^ Q j ' ^^^^ 

where we have defined the 2x2 matrix M. Eq. (fT2l) leads to 

Ao = MuA^N (13) 

with Mil denoting the element on the top left comer of the matrix M. Therefore, the relation ([13]) 
allows us to rewrite eq. (|6]) in the form: 



A* A, -|Miil2 ^'^^ 



Our results will be illustrated in Sec. HVlbelow. 



IV. MICROSCOPIC DISORDER AND THE EXTERNAL THERMOSTAT 

We characterize the presence of disorder in the sequence of barriers letting p{\)d\ be the 
probability distribution of the widths Aj of a generic barrier, with i = {1, A^}, to take values 
in a range d\ centered on A. In particular, we used pseudo-random generators to investigate two 
relevant choices for p(A). The first is the uniform density, with support on the unit interval, while 
the other is the density P = \ f e 2", supported on M+, which is obtained from the gaussian 
density p = Af{0, 1) by retaining only the positive values of the A/s . Each of the two distributions 



is characterized by the corresponding mean (A) and variance a [11911 . It proves useful to introduce. 
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for both these distributions, the realization mean and variance, defined, respectively, as 

N 

N 



i=l 
1 ^ 



i=l 

In the large limit, the random variable A^, which varies from realization to eralization fo the 
sequence of barriers, converges in probability to the mean (A), while the random variable 
converges with probability 1 to a^. Since we use dimensionless variables in eq. (|7]), we introduce 
the rescaled barrier width as: 

Xi = cXi, with c = — — . (15) 

il + (3)NXB 

Therefore, for any given and (3, the rescaled mean 

/3 



A/ 



(l + /3)iV 

is no longer a random variable, and attains the same constant value independently of the realiza- 
tion, hence on the density p. On the other hand, the rescaled realization variance 

2 ^ 
i=l 

remains a random variable which, for large converges to cr^ = c^a^ with probability 1. We 
introduce the vector- valued random variable An, defined by 

Aiv = {Ai,...,A7v} , (16) 

which corresponds to a given realization of the sequence of barriers and will be referred to 
as a microscopic configuration. For given (3 and A^, one may, then, consider the collection 
r2 = {Aj^\ Aj^''^} of Nr random realizations of the sequence of barriers which have been 
constructed numerically. Then, the average of a random observable O over the sample f2 simply 
reads as: 



Among the possible configurations, the regular one 

A^ = {A^...,A^} , (18) 



Nf = 1 00 - 
reg. lattice 




1 00 200 300 400 500 600 700 800 900 1 000 




900 1000 



FIG. 2: (Color online) Left panel: Behavior of J{N, Ajy, V, E; T) vs. N, for V = E = 9.60, Ey = 
leV, r = 10^, /3 = 0.1 and for different microscopic configurations: the red curve corresponds to a single 
uniformly distributed configuration Ajv, the blue curve represents the average {J)n over an ensemble of 
uniformly distributed configurations, whereas the dashed green curve corresponds to the regular lattice 
configuration ([TST l. Right panel: Same curves as those shown in the left panel, but here evaluated at E = 
5.16, where a condition close to the tunneling resonance is realized. 



m, will be 



which approximates the infinite superlattice of the literature on Bloch waves [|17 , 
crucial also in our work. 

In our numerical simulations we investigated the behavior of the coefficient J{N, A^, V, E; T) 
(we do not explicitly indicate the dependence on the parameters r, (3 and E^, not to overload the 
notation) at a given temperature T on the parameters of the model, in particular the number of 
barriers N, the height of the barrier V , and, mostly, the microscopic configuration A^v- 

Figure m illustrates the behavior of J{N,An, V, E; T), J{N, A^, V, E; T) and of the average 
{J)n, at T = 300K, over an ensemble characterized by p(A) = 1, for two different values of 
the energy E, one of which leads to a condition close to that of resonant tunneling. The plots 
show the behavior of J corresponding to a single disordered realization, of the average {J)n, 
computed over a set of A',. = 10^ disordered realizations, and of the value, hereafter denoted as 
Jb{N, V, E; T), corresponding to the regular lattice configurations (ITFi . They also reveal that (J)q 
tends, for growing N, to the most probable value of J{N, A^, V, E; T), given by Jb{N, V, E; T), 
as also discussed below. This holds for all values of V and of E, and was already observed 
in the equilibrium version of this model, discussed in nl6\ . Let us now consider a thermostat 
located at the left boundary, so that the plane waves entering the bulk have an energy obeying a 
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I ' ' 1 ' 1 ' 1 ' 1 1 I ' ' 1 ' 1 ' 1 ' 1 1 

200 400 600 800 1000 1200 1400 1600 1800 2000 200 400 600 800 1000 1200 1400 1600 1800 2000 

N N 

FIG. 3: (Color online) Left panel: Behavior of J{N, An, V; T) vs. for = 20, T = SOOK, = leV, 
r = 10^, corresponding to two disordered configurations drawn from tlie "iialf-normal" distribution (dark 
and light blue curves) and to the regular lattice configuration (thick green dashed curve). Right panel: Same 
as in the left panel, with two configurations obtained from the uniform probability density. 

classical equilibrium distribution at a given temperature T. In the following plots we consider a 
one-dimensional Maxwellian probability density 

and we average over all energies to obtain 

POO 

J(iV, A^, V;T)= / J(Ar, A^, V, E; T)f,,{E)dE (19) 
where the coefficient J, defined in eq. (fT4l) . is integrated in the r.h.s. of eq. (fT9] ). 

V. NONEQUILIBRIUM FLUCTUATIONS 

We show here some numerical results concerning the value of J for a single realization, aver- 
aged over the equilibrium distribution of energies, as shown in eq. (fT9l) . The typical behavior of 
J{N, Ajy, V; T) as a function of the number of barriers is illustrated in Fig. |3l the left panel 
shows the fluctuations, around the value Jb{N, V, E; T), of the transmission coefficient pertain- 
ing to two disordered configurations obtained from the "half-normal" distribution introduced in 
sec. |IVl Similarly, the right panel displays the fluctuating behavior of J with A^ for two random 
uniformly distributed configurations. 

Figured anticipates two further crucial aspects which will be addressed in more detail below. 
The first concerns the magnitude of the fluctuations of the values of J{N,An, V; T), Op{N, V) = 

10 
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FIG. 4: (Color online) Left panel: Behavior of J{N, An,V;T) vs. for = 20, T = 300K, = leV, 
r = 10^ in an ensemble of random uniform configurations for different values of V. Right panel: Behavior 
of J{N, An, V; T) vs. iV for F = 20, T = 300K, = leV , r = 10^ in an ensemble of random uniform 
configurations for different values of j3. 



a/ ((J — {J)n)'^)n, which decreases with A^. 

This decay of the size of the fluctuations allows us to identify a "mesoscopic" scale Nmeso- The 
number Nmeso depends on V and the coefficient J{N, An,V;T) depends only weakly on the 
microscopic configuration, if > Nmeso- A second scale N macro > N^eso is identified such that 
J{N, An,V]T) depends neither on the configuration nor on the number of barriers if > Nmacro- 
In the case A^ > Nmacro we practically have an infinitely finely structured sample of macroscopic 
fixed length. 

We also investigated the dependence of J{N, A^, V; T) on V and (3. The left panel of Fig. |4] 
corroborates, in the limit of large A^, the numerical results illustrated in Fig. |3] Namely, the trend 
of the random values J{N, A^, V; T) to approach Jb persists even when the potential energy of 
the barrier V , as well as the amount of insulating fraction in the system (related to /3) are changed. 
This stems as one of the prominent features of our model. In particular. Fig. S] shows that an 
increase in V ox'm (3 produces a decrease of the observed transmission coefficient. 

Next, in the left panel of Fig. [5l we plotted the behavior of {J)n{N, V; T) with A^ and we 
compared it with the value Jb{N, V; T) pertaining to a regular lattice. The result is consistent 
with those of Figs. [3]and|4l for it shows the regime corresponding to A^ ^ 1, where the curve of 
{■J)niN, V; T) varies very slowly with A^, approaching the value Jb{N, V; T). This convergence 
process is a collective effect, in that it is related to the de localization of the wave function described 
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FIG. 5: (Color online) Left panel: Behavior of {J)n vs. NforV = 20, T = SOOivT, = leV, /3 = 0.1 
and r = 10^, in an ensemble of random uniform configurations for different values of Nr. Right panel: 
Behavior of {J)n vs. N for Nr = 100, T = 300K, = leV, r = 10^ and for different values of V, in 
an ensemble of random uniform configurations. 



by the Bloch waves theory [|l8(] for regular lattices. It is also worth pointing out that Jb{N, V; T) 
resembles the value predicted by the Bloch theory only in the limit ^ 1, in which the model 
corresponds to a good approximation of the regular infinite lattice. The right panel of Fig. |5]shows 
the behavior of {J)q{N, V; T) for different values of V , and highlights the dependence of Nmacro 
on V. In particular, the figure shows that taking N > Nmacro ~ 150, 100, 50 for, respectively, 
V = 20, 60, 100, allows one to reach a good accuracy even with small samples f2. 
Moreover, we verified that the energies of the incoming wave, at T = 300K, lie in the conducting 
band of the infinite periodic chain of barriers, which implies that Jb{N, V; T) > 0. Furthermore, 
the plots in Figs. |6] and |7] show the behavior of J(N,An, V; T) for different values of and r. 
At fixed Ey, we see that varying r yields an increase of the transmission coefficient. Similarly, 
increasing E^ results, in the observed region of the parameter space, in a slight increase of the 
transmission coefficient. 

Figure[8l instead, shows the behavior of ( J)n, averaged over an ensemble of = 10^ random 
uniform configurations, for different values Ey. It is worthwhile noticing that, in the limit of 
vanishing external fields, the value of the transmission coefficient is sensibly close to the value 
reported in [|l6ll . referring to an equilibrium model. A few comments can be drawn, here, also 



based on the comparison with those results discussed in Ref. [|l6|] . In the first place, the absence 



of localization can be traced back to the fixed finite amount of insulating material, which we have 
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FIG. 6: (Color online) Behavior of J{N,An, V; T) vs. N with y = 20, T = SOO/C, = leV , p = 0.1 
and, respectively, r = 10 (top left), r = 10^ (top right), r = 10^ (bottom left) and r = 10^ (bottom right). 



even in the N ^ oo limit, because L is fixed. As a consequence, incoming waves may, at most, 
be damped by a finite factor, except, perhaps, for a negligible set of energies which we have not 
observed. This distinguishes our model from the tight-binding model, which is more extensively 
investigated in the specialized literature, and also prevents the application of the Furstenberg's 
theorem [11]. Indeed, introducing the (A^+ l)-th barrier in one of our system realizations produces 
a rearrangement of the previous barriers. Mathematically, this means that the product of the first 
random matrices is replaced by a new product. Differently, the case of ergodic-like theorems, 
such as Furstenberg's theorem, applies to products of A^ random matrices which do not change 
when they are multiplied by the (A^ + l)-th matrix. 

The decrease of the size of fluctuations with A^, which will be explored in more detailed below, can 
be regarded as a phenomenon of self-averaging of the observable J [|ll|l . In particular, our results, 
further supported by the analysis of the PDF of the transmission coefficient. Figs. l9l and [TOl below, 
show that the random values J(A^, A^v, V] T) converge in probability to Jb in the A^ — t- oo limit. 
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FIG. 7: (Color online) Behavior of J{N,An, V; T) vs. N with y = 20, T = SOO/C, = 2eV , p = 0.1 
and, respectively, r = 10 (top left), r = 10^ (top right), r = 10^ (bottom left) and r = 10^ (bottom right). 





FIG. 8: (Color onUne) Left panel: Behavior of (J)^ vs. E.^ for N = 300, F = 20, T = SOOJT, ;3 = 0.1 
and r = 10^, in an ensemble of Nj. = 10^ random uniform configurations. Right panel: Magnification of 
the framed part of the left panel. 
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As shown in Fig. [TOl given a sample 17 of uniform realizations, JsiN, V; T) corresponds to the 
most probable value of the random variable J in the sample, which, when grows, tends also to 
the mean {J)n{N, V; T). Let us now investigate, more accurately, the structure of the fluctuations, 
in the sample Vt of random uniform configurations at temperature T = 300/^. Denote by 
Pn{J) the probability density pertaining to the random value J. The numerical results presented 
so far on the relation between {J)q(N, V; T) and Jb{N, V; T), as well as on the decrease of the 
fluctuations size with growing A^, indicate that Pn{J) peaks more and more around the reference 
value Jb{N^ V; T). To show this more clearly, we numerically calculated the quantity Pn{J) and 
we plotted in Fig. |9]the resulting curves for different values of N and for V = 20. The maxima of 
the PDF in Fig. |9]are approximately located at J = Jb(^7 V', T), cf. also Fig. [lOl and tend, for 
large A^, to the mean value {J)q(N, V; T). One further realizes that Pn{J) obeys a sort of large 
deviation principle. By this we mean that the limit 

Cn{j) - ^ a-j) (20) 

exists for the collection of J values conceming the sample of different realizations. Figure [10] 
shows, for the range in which we have good statistics, that is apparently smooth and strictly 
convex like a normal large deviation functional. However, it is worth pointing out, again, that the 
N oo limit is not achieved in the standard fashion of products of random matrices. Moreover, 
even the observable J is not of the usual kind discussed in large deviation theory, in that it is not 
given by a sum of i.i.d. random variables, being it related, in general, to the random variables 
Ai, Aat in a highly nonlinear fashion. 

It is also interesting to note that the locus of the minima, i.e. of highest probability density, of 
the curves Cn{J) is represented by the locus of the values Jb{N, V; T), Fig. \T0\ Moreover, Fig. 
nn reveals that, for small deviations from Jb{N, V; T), Cn{J) is quadratic, as expected where the 
central limit theorem applies. 

From the validity of a principle of large deviations and of the central limit theorem, one expects 
the following asymptotic behavior for the fluctuations of the transmission coefficient: 

'-AT) , forAr»l. (21) 



We numerically checked that this is indeed the case by evaluating the ratio in eq. (|2T)) for V = 20, 
cf. the top panel of Fig. [121 The bottom panel of Fig. [I2l further illustrates the decrease of 
Nmeso with V: taking Nmeso such that crp{N, V) < a* = 5.0 x 10"^ for A^ > Nmeso, where a* 
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FIG. 9: Probability densities pat (J) for different values of N and for V = 20, T = 300 JT, = leV, 
/3 = 0.1 andr = 10^. 
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FIG. 10: Rate functional associated to pn{J) for different values of N and for y = 20, 1^ = 20, 
T = 300K, Ey = leV, /? = 0.1 and r = 10'^. The curves ^Af move downwards for growing N, so that, as 
expected, in the N ^ oo limit, ({J) intersects the horizontal axis only in J = J^. 
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FIG. 11: Rate functional Cn{J) associated to pn{J). Blue points: results of the numerical simulation with 
N = 100 and the same set of parameters used in Fig. |9] Green dashed line: Fitting of numerical data with 
the parabola a{J/JB - if + b, with parameters a = 1.25591 ± 0.2575 • lO^^ and b = -2.2858 • lO'^ ± 
5.36 • 10"^. 



is considered small, Nmeso{V) is found to rapidly decrease with V: we obtain Nmeso ~ 90, 55, 30 
for, respectively, V = 20, 60, 100. Then, Nmeso must tend to when the barrier height grows, 
because the transmission coefficient vanishes in this limit for any configuration A^v, hence the 
fluctuations are also suppressed. 



VI. CONCLUSIONS 



In this work we investigated a quantum mechanical model under nonequilibrium conditions, 
and focused on the role played by the disorder on the transmission coefficient J. Our numerical 
investigation reveals the existence of appropriate mesoscopic and macroscopic scales, respectively 
denoted by Nmeso and Nmacro, which are not as widely separated as in thermodynamic systems. 
The novelty of our approach stems, first, from the introduction of a thermal average of the trans- 
mission coefficient over an equilibrium distribution of energies guaranteed by the presence of an 
external thermostat at a given temperature T. Furthermore, we also proposed a novel approach 
to deal with the "thermodynamic limit" of the model: we prescribed a fixed, macroscopic, length 
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FIG. 12: Top panel: Fit of cTp/ {J)n (red curve), with the same set of parameters of Fig. |9] with the curve 
f{N) = ax N'K with log(a) = -5.13922 x 10~^ ± 4.605 x 10"^. The plot confirms the asymptotic 
behavior of Up/ {J)n, given by eq. (|2T]) . Bottom panel: Behavior of ap{N, V) vs. N, for V = {20, 60, 100}. 
The thick orange dashed line denotes the value a* below which we may conventially assume that the size 
of fluctuations is practically negligible. 

for the system, so that the increase of the number of barriers does not yield the divergence of the 
overall length, rather it results in a more and more refined partition of barriers and wells. The 
novel route proposed in this work leads, in general, to different results with respect to the standard 
thermodynamic limit discussed in the literature, and also makes some classical results, e.g. the 
Furstenberg's theorem, not applicable. Our approach might, hence, open a new line of investiga- 
tion in the theory of disordered systems and could also allow to shed new light on the transition 
from the microscopic to the macroscopic scales. Our numerical results suggest that, in presence 
of off-diagonal disorder, the wave function is delocalized over all the system length, thus no local- 
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ization effect, of the like typically occurring with systems perturbed with diagonal disorder \\l5n . 
occurs. Moreover, the disorder, at the microscopic level, induces an irregular behavior of the trans- 
mission coefficient J. The variable J is self-averaging for growing A^, and peaks over the most 
probable value of J. Interestingly, this value is Jb, which corresponds to a microscopic ordered 
array of barriers and conducting regions. Moreover, in the ^ 1 limit, large deviations from 
the value Jb are possible, and the structure of these fluctuations is governed by the rate functional 
C{J), which has been numerically determined. It is worth emphasizing that, at variance with the 
standard derivation of large deviation principles for the sum of independent and identically dis- 
tributed random variables, our results hold for a random variable J which is related, in general, 
in a highly nonlinear way to the random widths of the single barriers. Finally, this work suitably 
extends to a given nonequilibrium regime the results reported in [16], obtained in absence of the 
external field. 
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